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We report a cluster of results on fc-QSAT, the problem of quantum satisfiability for fc-qubit 
projectors which generalizes classical satisfiability with fc-bit clauses to the quantum setting. First 
we define the NP-complete problem of product satisfiability and give a geometrical criterion for 
deciding when a QSAT interaction graph is product satisfiable with positive probability. We show 
that the same criterion suffices to establish quantum satisfiability for all projectors. Second, we 
apply these results to the random graph ensemble with generic projectors and obtain improved 
lower bounds on the location of the SAT-unSAT transition. Third, we present numerical results on 
random, generic satisfiability which provide estimates for the location of the transition for k — 3 
and k = 4 and mild evidence for the existence of a phase which is satisfiable by entangled states 
alone. 
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I. INTRODUCTION 

The central achievement of classical complexity the- 
ory was the discovery that certain natural classes of con- 
straint satisfaction problems, prototypically boolean sat- 
isfiability (SAT), can encode all verifiable decision prob- 
lems in an efficient manner. This insight, encapsulated 
in the Cook-Levin theorem, and the general consensus 
that P 7^ NP, imply that any classical approach to solv- 
ing instances of SAT must in general take exponentially 
long pp. Unfortunately, the recent development of an 
analogous set of quantum complexity theoretic results 



suggest that a similar situation holds even for quantum 
computers: in this case, quantum satisfiability (QSAT) 
can encode any problem in the class QMAi [E] and ought 
therefore resist solution by even the most subtle quantum 
attacks [2]. 

These hardness results, while rigorous, are in some 
sense unsatisfactory. They imply that (quantum) satis- 
fiability is exponentially hard to solve in the worst-case, 
but shed little light on what separates the worst-case in- 
stance from the trivial. In this pursuit, it is useful to in- 
troduce a measure on instances of satisfiability and study 
the statistical features of the ensemble rather than the 
worst-case behavior of the whole class. This approach 
provides control parameters with which to tune our at- 
tention across various regimes of difficulty and also allows 
us to bring to bear the full arsenal of statistical mechanics 
and probability theory 02 14] . 

This paper is devoted to such a study in the context 
of quantum complexity theory. We study quantum satis- 
fiability (fc-QSAT), the decision problem of determining 
whether a given Hamiltonian H acting on N qubits has a 
ground state with precisely zero energy. The correspond- 
ing Hamiltonian has the form 

M 

ff=E r (!) 

where IT™ is a fc- local (i.e. fc-qubit) projector associated 
with the (hyper-)edge m = (mi, m%, • • • , rrik) of an inter- 
action graph G. Such a Hamiltonian has a ground state 
\ip) of zero energy if and only if \tp) is simultaneously an- 
nihilated by all of the H m . Furthermore, an instance of 
fc-QSAT must come with a promise 6 — l/poly(N) that 
the ground state energy of H, should it be greater than 
zero, will actually be greater than 8. We will have lit- 
tle more to say regarding the promise gap in this paper, 
although we believe that it is satisfied in the thermody- 
namic limit of the ensembles we consider. 

In Bravyi's original work introducing the fc-QSAT 
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problem, it was shown to be QMAi-complete[2] for k > 4. 
In a previous paper [5], we introduced a random ensem- 
ble for fc-QSAT and worked out the coarse features of 
its phase diagram in the large N limit. There are two 
sources of randomness in the /c-QSAT ensemble: 1. the 
discrete choice of interaction graph G and 2. the con- 
tinuous choice of projectors H m associated to each edge. 
The latter choice is particularly powerful: generic choices 
of projectors reduce quantum satisfiability to a graph, 
rather than Hamiltonian, property. More precisely, for 
fixed G, the dimension Rq = |kcriJ| of the satisfying 
subspace of H is almost always minimal with respect to 
the choice of projectors H m . Thus, if G can be frustrated 
by some choice of projectors, it is frustrated for almost 
all such choices. This "geometrization" property allows 
us to make strong statements about the quantum sat- 
isfiability of Hamiltonians associated with both random 
and non-random graphs and even non-generic choices of 
projector, both analytically and numerically. Through- 
out this paper, we use the term generic to refer to the 
continuous choice of projectors and random to refer to 
the choice of the graph. 

This paper presents a cluster of results regarding 
generic quantum satisfiability under the restriction that 
the projectors II™ 1 = \<f>) [4>\ arc rank 1. In Sec. [TT| we 
introduce and characterize product satisfiability (PROD- 
SAT) on non-random graphs G. This essentially classical 
property is NP-complete, but will allow us to prove that 
interaction graphs in which each clause (edge) may be 
matched with a neighboring qubit (node) are quantum 
satisfiable for any choice of projectors. For the marginal 
case M = N, in which these correspond to perfect match- 
ings of the clause-qubit graph, the associated product 
states are discrete and provide some access to the ground 
state space, as will be explored in Sec. |II A| As physicists, 
we refer to the qubit-clause matchings as dimer coverings, 
although computer scientists will likely prefer to think of 
perfect matchings and Hall's marriage condition. 

The remainder of the paper is especially concerned 
with mapping out the phase diagram of generic quantum 
satisfiability on large random interaction graphs. As in 
our previous work|5] , we focus on the Erdos-Renyi ensem- 
ble with clause density a, in which each of the potential 
(^) edges appears with probability p — aN/ ( fc ) . As one 
might expect, the ensemble enjoys a satisfiable (SAT) 
phase at low clause densities a and becomes frustrated 
(unSAT) at sufficiently high clause density. Originally, 
we showed the existence of the SAT phase by construct- 
ing satisfying product states using something akin to a 
transfer matrix approach. This algorithm works generi- 
cally so long as the graph G does not have a hypercore - 
that is, a maximal subgraph in which every qubit partici- 
pates in at least two clauses. This property is guaranteed 
for a < aih c (k), the critical density beyond which a hy- 
percore appears. These critical values are available in 
the literature [6]; for k = 3, a.hc ~ 0.81, while for k — >• oo, 
a hc -> 0. 



state satisfiability of Sec. [TT]to the random graph ensem- 
ble. This improves the lower bound on the SAT/unSAT 
transition from to a^c, the critical value for the ex- 
istence of dimer coverings. For k — 3, ctdc ~ 0.92 is not 
much of an improvement over ahc ~ 0.81, but as k — > oo, 
ctdc — > 1 instead of 0. Moreover, this is a complete char- 
acterization of the phase in which the ensemble is generi- 
cally satisfiable by product states. If there is a satisfiable 
phase for larger a, its satisfying states will generically be 
entangled. 

At high clause densities, a simple local bound on the 
generic dimension of the zero energy subspace provided 
the first evidence for the existence of an unSAT phase. 
More recently, Bravyi et al (BMR) [7] have dramatically 
improved this bound by finding the ground state de- 
generacy of larger clusters ( "sunflowers" and "nosegay" 
graphs) and determining their prevalence within the large 
random graph. Using these techniques they have placed 
an upper bound, a+, on the SAT/unSAT transition of 
aj" ~ 3.594 for k = 3 and with an exponential scaling 
a+ ~ 2 fc for large k. 

The gap between the best rigorous lower and upper 
bounds on the SAT/unSAT transition is distressingly 
large: a~ = oy c = 0(1) versus a+ = 0(2 k ). Moreover, 
the existence of a satisfiable regime above o^e would pro- 
vide a large ensemble of large graphs with intrinsically 



In Sec. Ill we apply the characterization of product 



entangled satisfying states. Thus, in Sec. IV we address 
the situation for larger values of a by adducing some pri- 
marily numerical evidence about the existence of random 
instances which are SAT but not PRODSAT. We note 
that recently, Ambainis et al [8] have proven the exis- 
tence of such an entangled SAT phase for k > 12 using 
their newly proven quantum Lovasz local lemma coupled 
with the characterizations presented in this paper. We 
will have more to say about this exciting development in 
the conclusion. 



II. PRODUCT SATISFIABILITY ON GENERAL 
GRAPHS 

One of the remarkable results in Bravyi's original work 
on QSAT was the explicit construction of satisfying prod- 
uct states for k = 2, thanks to which 2-QSAT turned out 
to be easy to decide and the concomitant satisfying states 
easy to construct. This poses the question whether there 
exist regimes for k > 3-QSAT in which analogous results 
hold. We therefore define: 

Definition 1. An instance of QSAT, (G,n m ), is 
PRODSAT if it has a satisfying product state. That is, 
3 |*) = \ip x ) <g) • • • ® \ipff) s.t. IT™ I*) = Vm e G. 

PRODSAT may be viewed as a decision problem in its 
own right and, in the presence of a polynomially small 
promise gap S, it is efficiently verifiable by classical com- 
putation. A witness is simply a collection of 27V C com- 
ponents of a satisfying product state, the energy of which 
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FIG. 1: (Color online) (a) Example of a k — 3 interaction 
graph with M < N. Circles (green) indicate qubits and 
squares (red) indicate clause projectors that act on adjacent 
qubits; (b) a dimer (shaded blue) covering that covers all 
clauses. 

may be evaluated in linear time. Moreover, by choos- 
ing II m to project onto computational basis states, it is 
clear that PRODSAT contains SAT and is therefore NP- 
complete. 

In this section, we investigate product satisfiability for 
fixed graphs G and general choices of rank 1 projectors 
IF" = \4> m ) (4> m \. We show that G supports satisfying 
product states on open neighborhoods in the projector 
manifold if it has a qubit-clause dimer covering that cov- 
ers all of the clauses - ie. an assignment of a unique qubit 
to each clause (see Fig. [T]) . Appealing to the algebraic 
properties of the space of projectors, this shows that G is 
in fact PRODSAT for all choices of IP™. It also suggests a 
counting criterion for the ground state degeneracy. Thus, 
we will show 

Theorem 1. If an interaction graph G, viewed as a bi- 
partite factor graph of qubits and clauses, has a dimer 
covering of its clauses, then 

1. there exist open neighborhoods U C (CP 2 _1 ) M of 
the rank 1 projector manifold such that IP™ £ U => 
(G,FF") is PRODSAT. 

2. Moreover, (G,IF n ) is PRODSAT for all IF" e 

As a practical matter, product satisfiability is a ques- 
tion about the existence of solutions to a collection of 
M homogeneous equations in N CP 1 degrees of freedom 
(unless otherwise noted, paired i indices contract, but m 
indices do not): 

C-., t «-C t =0 VmeG (2) 

where ipm 1 is the ii'th component of the state on the 
mi'th qubit and 4> m * is the state onto which FF" projects. 
Since we are interested in generic choices of the 4> m , 
one might expect naive constraint counting to determine 
whether (M < N) or not (M > N) a given graph G 
is generically PRODSAT. Unfortunately, things are not 
quite so simple because a G with a low average density 
(M < N) may still have subgraphs that are dense enough 
to overconstrain the equations. 



To make this last observation sharper, let us briefly 
restrict our attention to product projectors. In this case, 
Eq. ([2| reduces to 

(C'ftXCt) • • • (C*^mJ = V m e G (3) 

This equation is clearly satisfied if any of its factors is 
satisfied, which requires at least one of the k qubits to be 
(uniquely) fixed orthogonal to the local projector. Thus, 
for generic <j> of product form, in which the local projec- 
tors on a site n are not parallel, the full set of M equa- 
tions is solvable if and only if one can uniquely associate 
a qubit n to each clause m. Graphically, this is equiva- 
lent to the existence of a dimer covering of the bipartite 
factor graph of G which covers all M clauses. Note that 
we refer to matchings as dimer coverings if each clause 
is uniquely associated to a qubit even if some qubits are 
left unpaired. See Fig [I] 

We now relax the constraint that <f> take product form 
and argue that the existence of the above dimer covering 
is still the relevant constraint counting for product satis- 
fiability. Fet us consider an instance (G, 4>) and assume 
that it is PRODSAT; w.l.o.g. we may choose local bases 
such that |\&) = (i) ® • • • ® (?) is a satisfying assignment. 
Choosing stereographic coordinates ( z ™) for each qubit, 
we find the constraint equation is: 

6 m - ■ z 11 z 12 ■ ■■ z ik = (4) 

where now superscript i\ indicates exponentiation and 
the i indices run from to 1. Since z n — satisfies this 
equation, we find that 4>qq...q = necessarily but that the 
remaining components of the 4> m are unconstrained. 

Now we perturb </> m i-> cj) m +8(f> m and attempt to follow 
the solution z n = i— > z n + 5z n . To linear order, we have 
(for each clause m E G): 

(5) 

This is a collection of M sparse linear inhomogeneous 
equations for the N variables 5z n . Fet Vt mn be the MxN 
matrix of coefficients on the left hand side. This matrix 
connects qubits to clauses and is sparse in the same pat- 
tern as the node-edge adjacency matrix A mn of G. 

Equation ([5| is solvable for arbitrary (5</>^ if and only 
if the matrix Q is surjective. Trivially, this requires M < 
N, but it also requires that some M x M subdeterminant 
of n be nonzero. Fet us define a polynomial discriminant 
of surjectivity as the tuple of all possible subdeterminants 
offi: 

disc r (fi) = ^ (-l)°"Oi T ( g (i)) • • • Omt(<t(m)) (6) 

where r : M '—t N runs over all injective maps. For for- 
mal components <f), disc T (!7) is not identically zero if and 
only if one of the terms in the above sums is not identi- 
cally zero; each non-zero term corresponds to a particu- 
lar pairing of clauses and qubits in the adjacency graph 
which covers all M clauses but need not cover all of the 
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nodes. Thus if G has such a dimer covering, Q is almost 
always surjective, whereas it is never surjective otherwise. 

We can summarize the above analysis as follows: for 
G having dimer coverings, generic choices of 4> for which 
(G, <fi) is PRODSAT may be extended to open neighbor- 
hoods around (j> which are PRODSAT. Moreover, G has 
product states at generic (non-parallel) product projec- 
tors. One can repeat the perturbative analysis of Eq. ^ 
near such product projectors and show that the dimer 
covering condition allows the same extension onto open 
neighborhoods in the full projector space. This proves 
part 1 of theorem [I] [TB] 

The above analysis is essentially local, showing that the 

set W C (CP 2 * -1 )*' of PRODSAT projector choices is 
nonempty and full dimension; to extend product satisfia- 
bility to the entire projector manifold, we need to appeal 
to some basic results in complex projective geometry: 
1. Any Zariski-closed full dimension subspace of an ir- 
reducible complex projective space (such as (CP 2 ~ )* f ) 
is in fact the whole space and 2. W is Zariski-closed. 
This shows part 2 of theorem [T] For more details, see 
Appendix [A] 

Finally, we note that the converse statement also holds 
(also shown in Appendix [A]) : 

Theorem 2. If G does not have dimer coverings, 
(G, IT™) is not PRODSAT for almost all IT™. 



they will remain so on an open neighborhood of 4>q, since 
the determinant of their overlap matrix is a smooth func- 
tion. We recall that the QSAT dimension Rq takes its 
minimal value almost everywhere on projector space [5]. 
Since Rps lower bounds Rq on an open neighborhood, 
it must therefore lower bound Rq for all projectors </>. 

Using the dimer covering states to characterize the full 
SAT subspace ker(iJ) for a given graph G requires an- 
swering two questions: 1. Are the product states linearly 
independent? and 2. Do they span ker(H)? For some 
simple families of graphs, it is possible to prove the lin- 
ear independence of the dimer states by exploiting the 
geometry of their dimer coverings. On the other hand, 
it is easy to construct graphs where the dimer states are 
not independent: for a fully connected bipartite factor 
graph, the number of dimer coverings is Nl, which is sig- 
nificantly greater than the size 2 N of the Hilbert space. 

The second question is harder to address analytically. 
Product states do span the kernel for many graphs that 
we have studied numerically but we also know examples 
of small graphs for which they do not. Rather, Rq can 
be strictly greater than the number of dimer coverings. 
Indeed, H restricted to product projectors may have an 
even greater satisfying dimension than it does for general 
entangled projectors, even though the number of product 
states should be no greater. 



A. Counting product states at M = N 

We now comment on the marginal case when M = N. 
In this case, one expects the satisfying product states, 
when they exist, to be discrete. Let us return to projec- 
tors of product form and consider a choice <po in which 
no local projectors on a single site n are parallel, as men- 
tioned above. In this case, it is clear that there is a 
one-to-one correspondence between dimer coverings c of 
G and product state solutions ip c . The number of prod- 
uct states at the special point <po is given by the number 
of dimer coverings. 

Now, one can pick any of the ip c and repeat the lo- 
cal perturbative analysis around cf>Q. This will show a 
unique way to continuously extend ip c in an open neigh- 
borhood of 4>o. Thus, at least on open neighborhoods 
within projector space, the number of product states is 
bounded below by the number of dimer coverings. Since 
the problematic points in c/> space for this local analysis 
are actually degenerate points with higher satisfiability, 
this should be a lower bound on the number of prod- 
uct states for any </> and it should be the exact counting 
almost everywhere. 

With a handle on all of these product states, one might 
hope to make some sharp statements about the QSAT di- 
mension of G - that is, the generic Rq — \ kei(H)\ for the 
graph G. Let PS — Span{i/> c } be the span of the dimer 
covering states at some generic product projector point 
4>q. If Rps = \PS\ of the ip c are linearly independent, 



III. APPLICATION TO RANDOM GRAPHS 

From the previous section, we know that the existence 
of dimer coverings on G implies quantum satisfiability 
- and that their nonexistence shows satisfying product 
states are no longer generic. For the random graph en- 
semble, we can therefore lower bound the QSAT transi- 
tion a c by the threshold for the presence of dimer cov- 
erings adc- This turns out to be precisely the point at 
which the hypercore of the random graph reaches the 
critical density Mh c = N) lCl beyond which it is certainly 
impossible to dimer cover the clauses. In general, the 
hypercore emerges at ay ic < 1 with a clause density less 
than I and as a grows, it gets denser. Only at adc > cthc 
does it become critically constrained. See Table [T] for 
numeric values of adc- 

In order to show that dimer coverings exist below 
adc, we rely on a somewhat indirect argument based on 
the results of Mezard et aL[9] regarding the random k- 
XORSAT problem, or, in physics parlance, the p-spin 
Ising glass (with p = k). In this problem, N classical 
bits live on an interaction graph G whose (hyper-)edges 
define parity check constraints. That is, for each edge to, 
we must choose J m G {0,1} and then determine whether 
the following linear constraint equations are satisfiable: 

■^mi ~t~ %rri2 ~t~ " " " ~t~ •^m^ — J ra (7) 

where the bits and all arithmetic takes place in 

the field Z2. 
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PRODSAT UNSAT (k = 3) 

0.17 0.81 0.02 3.59 1.26 7.49 

FIG. 2: Phase diagram of random fc-QSAT for k = 3. The 
labeled transitions are (left to right): the emergence of the 
giant component (oe gc ); the emergence of the hypercore (aha)', 
the disappearance of product states (a ps = Ode)', the best 
current upper bound on satisfiability due to Bravyi et al 7 
{at)\ the classical cavity transition for fc-SAT (a cav ) and the 
weak Pauling bound (a w b). At large fc, the known upper 
bounds all scale exponentially while the best lower bound 
(a ps ) is bounded above by a = 1. 



This defines a linear system of inhomogeneous equa- 
tions of exactly the same form as what we discovered 
when linearizing PRODSAT: 
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TABLE I: Summary of critical values for various random k- 
(Q)SAT properties. The values for the emergence of a hyper- 
core othc and for the dimer coverability/product satisfiability 
cip S = ot.de of the random interaction graph may be found in 
[9]. otj (3) comes from the BMR 'nosegay' bound[7] while the 
values at larger k and in the large k limit come from the 'sun- 
flower' bound (cf. Appendix [cj . For comparison, we include 
the best known rigorous bounds on the classical satisfibility 
transition [101111] . 
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where A is the node-edge adjacency matrix of G. In 
this case, A is surjective if and only if all choices of J m 
result in satisfiable instances of XORSAT. Moreover, if 
A is surjective over Z2, so is the sparse matrix £l mn = 
4> m nA mn over C with arbitrary components </> mn : if f2 
were not generically surjective, the rank discriminating 
polynomials disc T (f2) of Eq. [6] would be identically zero, 
including when mapped homomorphically to Z2. 

Mezard et al show that the random XORSAT problem 
is unfrustrated for a < a>d c - Indeed, they show precisely 
that with probability going to 1 in the thermodynamic 
limit, the random graph G is XORSAT for all choices of 
J m . In particular, this implies that A is surjective for 
a < ctdc and therefore that fl is as well, dimer coverings 
exist and the QSAT problem is satisfiable. Meanwhile, 
for a > add M > N on the hypercore and the dimer 
coverings are gone. 

We note that the previously known lower bound was 
given by a^c, the threshold for the emergence of a hyper- 
core in G. As fc becomes large, a^ c approaches zero while 
otdc approaches 1 (from below). This significantly raises 
the known lower bound for quantum satisfiability, but it 
still does not approach the best known upper bounds, 
which scale with 2 k . See the phase diagram in Fig. [2] for 
the numerical values at fc = 3. 



IV. BEYOND PRODSAT 

In the previous section, we argued that random fc- 
QSAT is PRODSAT for all k, at clause densities a < 
adc(k) = a ps (k) with a ps (2) = |, a ps (3) = .92... and 
a ps (k — > 00) — > 1~. (A heuristic way of understanding 
the latter limit is provided by the observation that for 
large fc, the probability of having sites with low coordi- 
nation is suppressed exponentially ~ e~ ak , so that the 
density of the random graph and its 2-core are almost 
the same.) 



What happens for a > a ps (k)7 Reliable information 
on this only exists for fc = 2: as the ensemble goes 
unPRODSAT, it immediately goes unSAT as well. For 
fc > 3, no such result has as yet been established. Rather, 
the best upper bounds on the extent of the SAT phase, 
due to BMR[7j, are well separated from the PRODSAT 
transition: a detailed analysis of the equations they de- 
rive yields the upper bounds of Table [I] The large-fc 
asymptote of this bound is a+(fc) = 2 fc_1 In 2 + O (jjx), 
which is exponential in fc, leaving plenty of space above 



1- for a SAT-unPRODSAT phase. 



In the following, we adduce some numerical pointers re- 
garding the existence of such a phase, which is arguably 
among the most interesting open questions in this prob- 
lem. 



A. Direct search for SAT— unPRODSAT instances 

First, we have directly searched for instances - irre- 
spective of their probability in the random graph ensem- 
ble - which are SAT-unPRODSAT. Our search for sat- 
isfying assignments proceeded via (i) exact diagonalisa- 
tion of systems with up to 14 qubits and (ii) a numerical 
nonlinear equation solver for finding product solutions. 
Our results were all consistent with the picture of prod- 
uct states outlined above. In addition, we found graphs 
with satisfying assignments with M > N . In these cases 
method (ii) failed to turn up any product ground states. 
We also diagnosed directly the presence of non-product 
ground states by considering the rank of the reduced den- 
sity matrix corresponding to a ground state after an arbi- 
trary portion of the system is traced over, cf. Fig. [3] The 
resulting rank of the density matrix is bounded above 
by the number of zero-energy states, provided these are 
spanned by product states (cf. Appendix|B for a proof). 
As an illustration, we show the corresponding rank for 
the 2-SAT problem in Fig. [3^a) for reduced density ma- 
trices of up to five qubits from instances containing 13 
qubits. While the SAT instances show the expected rank 
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a) k=2 



b) k=3 



k=3 



k=4 




1 2 3 4 5 
b qubits in block 



4 8 12 16 20 24 28 32 

Rank of 5 qubit RDM 



FIG. 3: (Color online) Rank analysis of reduced density matri- 
ces (RDM) in a ground state: (a) k = 2 : Rank of the reduced 
density matrix containing b — 1,...,5 qubits for different 
types of interaction graphs, (b) k = 3 : Histogram of ranks of 
6 = 5 reduced density matrices for all 252 different partitions 
(N = M — 10) for three different instances. 



for product states, the unSAT instance has full rank, sug- 
gesting generically entangled wavef unctions in the k = 2 
unSAT phase. In Fig. (3^b) we performed a similar analy- 
sis for k = 3. We determined rank histograms for all 252 
five-qubit reduced density matrices of N — M — 10 sys- 
tems for the three different instances listed in [17] . The 
PRODSAT instance leads to a histogram where the rank 
is bounded by the number of satisfying ground states 
(16 for this instance). The SAT-unPRODSAT instance 
shows ranks which are both smaller and larger than the 
number of satisfying ground states (23 for this instance), 
but do not reach the maximum of 2 5 = 32. Finally, the 
unSAT instance has full rank 32 on the two-core (the 
rank 8 occurring for some partitions is due to a satisfied 
dangling clause). We do indeed find that, in the presence 
of M > N on a (sub)graph, entangled states necessarily 
appear in any basis of the ground states. 



B. Finite size scaling 

Where is the the SAT/unSAT transition a c ? Numeri- 
cal studies on finite-site systems offer limited information 
as the properties of small random graphs are quite differ- 
ent from those of large ones. We do consistently find the 
presence of SAT instances with N > M: for 3-SAT, we 
investigated 101 random graphs for sizes with up to 13 
sites using complete diagonalization and up to 20 sites us- 
ing an iterative diagonalization technique, see left panel 
of Fig. [3] While complete diagonalization is limited by 
the growing matrix size but otherwise numerically stable, 
the iterative Lanczos procedure faces challenges when re- 
solving small energy gaps occurring at the SAT-unSAT 
transition. This leads to some instances becoming unde- 
cidable, i.e. not converging within several thousand iter- 
ations. The presence of undecidable instances is denoted 
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FIG. 4: (Color online) Probability to find SAT instances 
among 101 random graphs for (a) k = 3 and (b) k = 4. 
Numerical results obtained by complete diagonalization for 
N < 13 and iterative diagonalization for N = 14, 16 (k = 3, 4) 
and 18,20 (k — 4). For the latter system sizes filled sym- 
bols indicate the occurrence of (numerically) undecidable in- 
stances. 
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FIG. 5: (Color online) Finite size scaling of a c (N) for k = 
3, 4 obtained from the data displayed in Fig. [4] The empty 
symbols denote a values where the SAT ratio is 0.5. The 
lower (upper) ticks correspond to the largest a values with 
SAT ratio > 0.9 (smallest a values with SAT ratio < 0.1). 



by filled symbols [TH]. For 3-SAT, we find that curves 
for larger N nicely converge, suggesting a c sa 1 ± 0.06, 
where the error bar is estimated from the largest N in 
Fig. [5j For 4-SAT, instances are SAT for larger values 
of A = M — N (for fixed N), even though A appears 
to decrease with N (inset in right panel of Fig. [4 ) . The 
finite size scaling shown in Fig. [5] is also compatib e with 
a c w 1, but with a larger uncertainty of w 0.2. 

Thus, we find many examples of generically SAT- 
unPRODSAT graphs G and even some evidence that the 
unSAT phase transition in the random ensemble lies at 
a > a vs . 



V. IN CLOSING 

In this paper we have reported considerable progress in 
understanding product satisfiability. Applied to the ran- 
dom graph ensembles introduced in [S] , it yields a defini- 



7 



tive lower limit to the clause density at which satisfiabil- 
ity necessarily requires entangled states and certainly an 
improved lower bound to the location of the SAT-unSAT 
phase transition. Between the work reported in this pa- 
per and by BMR [7], there has thus been a fair amount 
of progress in narrowing the location of this phase tran- 
sition from the initial work [S]. 

Nonetheless, there remains a large gap between the 
best lower bound a~ — a ps = O(l) and the best up- 
per bound a+ = 0(2 k ). By contrast, the classical ran- 
dom /c-SAT problem has tight bounds on its satisfiability 
transition that scale as 2 k [10]. It is not obvious how to 
translate the techniques used in proving the tight lower 
bounds for the classical problem to the quantum case. 

Ambainis and co-workers [8] had recently proven a 
quantum version of the Lovasz local lemma, which di- 
rectly shows the quantum satisfiability of graphs G of 
bounded degree. Due to this degree restriction, the 
lemma does not directly apply to the Erdos random 
graph ensemble. Thus, they have used the characteri- 
zation of product satisfiability shown here coupled with 
the quantum Lovasz lemma in order to show that a wide 
SAT-unPRODSAT phase in the parameter a exists: the 
SAT regime extends to a c > 2 fc /(12 • e • k 2 ) which, for 
k > 13, is significantly beyond the PRODSAT transition 
at a ps ~ 1. 

Many other questions remain open and it is to these 
that we turn in closing: 

1. Does the SAT-unPRODSAT phase exist for all k > 
3 or does it really only arise at some larger fc? 

2. Does an easily checked classical graph property 
encode generic quantum satisfiability? Does a 
straightforward function compute the generic di- 
mension Rq7 

3. What can be said about higher rank projec- 
tors? The product satisfiability condition trans- 
lates straightforwardly to the higher rank case by 
simply treating a rank r projector as a stack of rank 
1 projectors. Thus, higher rank QSAT instances 
may naturally yield more entangled ground states 
than the rank 1 case. Indeed, preliminary work 
indicates that this is the case|12j. 



Acknowledgements 

We thank S. Bravyi, C. Moore, T. Osborne, D. Nagaj, 
A. Obus and especially O. Sattath for many enlightening 
conversations. We would also like to thank the generous 
hospitality of the Erwin Schrodinger Institute in Vienna 
where parts of this work were completed. 



Appendix A: Formal PRODSAT proof 

We expand on the complex algebraic details of theo- 
rems [l]and[2] That is, we show that interaction graphs 
G with dimer coverings are PRODSAT for any choice 
of projectors {II m }^f =1 and conversely that G without 
dimer coverings are not PRODSAT for almost all such 
choices. The relevant algebraic geometry background can 
be found in tUm]. 

Proof. Let P = (CP 2 ^ 1 ) M be the space of fc-qubit pro- 
jectors, S = (CP 1 )^ be the space of product states, and 
V C PxSbe the space of pairs of ({II m }^ =1 , {|V>n)}£Li) 
such that 1^) = \ipi) ® • • • Cg) IV'jv) satisfies the projectors 
II m . We note that P, S and P x S are each irreducible 
complex projective spaces in the Zariski topology and 
that V, as the zero-locus of collection of homogeneous 
algebraic constraints, is a closed subspace of P X S. Let 
7r : P x S — > P be the projection map onto the first 
component. Since S is projective, this map is closed. 

We are interested in the set of projectors W C P which 
are satisfiable by product states. These are precisely the 
projectors in the image of V under the projection tt: 

W = 7r(F) (Al) 

Since it is closed, W is closed. By the local arguments 
of Section [IlJ we know that W is full dimension and 
nonempty. Since P is irreducible and projective, it has 
no proper full dimension closed subspaces; ie. W = P. 

Conversely, if G does not have dimer coverings, there 
exist choices of product projectors H m which are clearly 
unsatisfiable by product states (see Eq. (|3|). Hence W 
is a closed proper subspace of P and cannot be full di- 
mension. Thus, W has codimension at least 1 in P and 
almost all choices of projector H m will fail to be product 
satisfiability. □ 

Appendix B: Reduced density matrix diagnostic for 
product states 

Let us consider a graph A and a wavefunction 
which is expressible as a sum of wavefunctions each of 
which is a product wavefunction on the sites i e A. Let 
7 C A, and let p = |*)(*| and p 7 = Tr A \ 7 p. Then 

rank p 1 < < Nq where No is the number of linearly 
independent product states, and Nq 1 ^ the dimension of 
their span when restricted to 7. 

Let us label a product basis for the ground state 
manifold by: {$ a \a = 1 . . . N } with |$ Q > = (gfe^?), 
(<f>?\4>?) = 1- Here ||A|| = n s . Thus, ($ a |*o) = 1 (nor- 
malised) and = nr=i (&? l&i ) ( not necessarily 
orthogonal) . 

A general ground state is 

1 n 

i*) = s E A *> ( B1 ) 

a=l 



with 



N n s 
0,5=1 i=l 



Tracing over the qubits on subset A \ 7 of the set A 
yields: 



(B2) 



for which 



j N 



(B3) 



JV 



5^ A^A Q [^ £7 |^) 



a,P=l 



n 

fc'eA\7 



/5 7 is at most an No x A/q matrix, with at most rank Nq irrespective of the choice of 7. 

If dim (span {(g)j' g7 |()f)")|a = 1 ... A^o}) = A^ 7 ' < Ao, then we can write, for a — Nq (say): 
Z)o=i 1 ^" V ° ) [®*'6'rl^?)]) so tnat terms like 



N -l 



E ^v; (JVo) [®^i^>] 



(B4) 



(B5) 



take on the same form as the ones in (B4). Using such 



substitutions repeatedly, one finds that the reduced di- 

(7) 

mensionality manifestly has a rank of at most Nq . 

In particular, for the case of 2-QSAT and A being a 
tree (||A|| + 1 = Nq), N^ = || 7 || + 1. When a qubit 
being traced out is a dangling one, this result is obvious. 
In the case of an internal qubit (labelled 2, say), we note 
that the Bravyi construction on A induces one on 7 via 

— T12T23, where represents the Bravyi transfer 
matrix between for the projector ij, so that the resulting 
state counting is that of the smaller graph. 



Appendix C: Improved estimate of from BMR 
sunflowers 



In this section we provide an integral representation of 
the upper bound on the ground state dimension obtained 
in [7] by studying 'sunflower' graphs for generic k. By 
means of this we can obtain upper bounds a+(k) for k > 
4 and an asymptotic estimate for large k. From [7] we can 
read an upper bound to the dimension of the satisfying 
space Rg as 



d=0 

where 



lim ±R G <lii2 + f2a d (d ln(l - 2~ k+1 ) + In (l + ^^)) = S, (CI) 



a,, = j 1 dt {kat " ^ e-^'' = t dt a d (t). (C2) 



d\ jo 



By inverting the sum over d and the integral over t we find 



S = In 2 + aln(l - 2~ fe+1 ) + jf dt(\n (l 



2 k - 2 / 

The last is a ^-dependent average with the Poisson distribution adit). We use the representation 



00 e - s 



(C3) 



ln(l + x)= / ds (l-e- sx ) (C4) 

s 
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to perform the average over d first (using the generating function of the Poisson distribution) and the integral over t 
afterwards. We are left with 

5 = i» 2 + t ,Mi-2-*«) + /^s?f(i-^(r( f l I )-r( J i T ,,))) (C5) 

where q ~ ka(l — e^ 5 ^ 2 *^ 2 )) and T(z,x) is the incomplete gamma function (with r(z,0) = T(z)). This integral 

representation, although seemingly complicated, is perfectly suited for numerical integration, since the integrand is 

always finite and cut-off at s = O(l), beyond which it decreases exponentially in s. 

I 



One can obtain the numbers quoted in the text for the exact in the large k limit. 



threshold a^(k) to arbitrary accuracy using Mathemat- 
ica. Moreover in the limit k, a — > oo and a/2 k finite, 
relevant for the large k asymptotics, we can bound the 
value of the integral, solve for S = and find the quoted 
result 



= 2 



fel ^ + 0(l/lnfc). 



(C6) 



it is worth noting that this bound is half of the value 
of the lower bound of fc-SAT which is known to become 



We have not been able to analyze in the same detail 
the upper bound given by the 'nosegay' graphs of [7J. 
Although this bound would be stricter for any k than 
that coming from sunflower graphs, our preliminary re- 
sults point towards the same asymptotic scaling ( C6 ) . It 



would be interesting to understand whether for large k 
this bound is tight, as happens in classical random fc- 
SAT. 
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